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Boundary effects in liquid-crystalline phases can be large due to long-ranged orientational correla- 
tions. We show that the chiral hexatic phase can be locked into an apparent three-dimensional N+6 
phase via such effects. Simple numerical estimates suggest that the recently discovered "polymer 
hexatic" may actually be this locked phase. 



PACS numbers: 64.70.Rh, 61.25.Hq, 87.16.Ka 

Liquid crystals provide a tabletop laboratory to study 
broken symmetries and the resulting low-energy Nambu- 
Goldstone modes Viewing two-dimensional crystal- 
lization in this way led to the proposal of an intervening 
phase between the liquid and the solid: the hexatic [0 , a 
two-dimensional phase with long-range orientational or- 
der but short-range positional order. In three dimen- 
sions, arguments based on Landau theory suggest that 
there is no similar orientationally ordered phase between 
the three-dimensional liquid and solid ||. However, be- 
tween a nematic phase and a crystal phase, the mesogens 
can have hexatic order in the plane perpendicular to the 
nematic director [Q. While this "N+6" phase could, in 
principle, exist (it is not terribly different from a biaxial- 
nematic phase) it was not immediately discovered. Re- 
cently, however, Strey, et al. Q have examined a phase 
of DNA, having a structure consistent with N+6 order, 
which they call the polymer hexatic. In this letter, we 
will consider the effects of boundaries on an N+6 phase 
consisting of chiral molecules, such as DNA. We will ar- 
gue that while symmetry predicts that the hexatic bond 
order should twist || surface effects can "lock-in" a pre- 
ferred set of directions which can effectively unwind the 
twisting hexatic. 

The phase sequence of a chiral mesogen usually starts 
with the cholesteric at the highest temperatures or low- 
est concentrations. As the temperature is lowered or the 
density increased, the mesogens can form any number 
of phases, including the (achiral) smectic-A phase, the 
smectic twist-grain-boundary (TGB) phase Q or the chi- 
ral hexatic ||. In the last case, the free energy which 
governs the transition from cholesteric to chiral hexatic 
depends on both the nematic director n and the complex 
hexatic order parameter tpQ-. 

F hu i k = J d 3 x j|(V - tq Q n) ^ 6 | 2 +r|^ 6 | 2 +u|V> 6 | 4 

+-K 1 (V ■ n) 2 + -K 2 (n • V x n + q f 

+±K 3 [n x (V x n)] 2 j (1) 

where r and u are Landau parameters and Ki are the 
Frank elastic constants. In the cholesteric phase, ( ipa ) = 



and (P reduces to the free energy of a cholesteric 
with equilibrium pitch P = 27r/q - When hexatic or- 
der persists yb & = |?/; 6 |e 6 and, as in the Meissner 
phase of superconductors the nematic director sat- 
isfies V x n = 0. However, in this phase, it is straight- 
forward to see that the hexatic bond-order rotates about 
the average nematic direction: V$6 — <Zo n o- Recalling 
the phenomenology of the superconductor or smectic-A 
liquid crystal J^] we emphasize that if the director un- 
winds then the hexatic bond-order must twist. Further, 
we expect the cholesteric pitch to unwind continuously to 
infinity at the cholesteric-to-chiral-hexatic phase transi- 
tion, as has, in fact, been observed ||. Possible phase 
diagrams are shown in Figure El 
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FIG. 1. Phase diagram of a chiral mesogen that can form 
an N+6 phase. In the cholesteric phase there is no hex- 
atic order (ip6 ~ 0) and V x n 7^ 0. In the chiral-hexatic 
phase V x n = but ^6 7^ and V#e = <?on. (a) Type-I 
behavior with no intervening defect phase, (b) Type-II be- 
havior in which a twist-grain-boundary phase (akin to the 
Abrikosov flux-phase in a superconductor) intervenes between 
the cholesteric and chiral-hexatic phase. 

Deep in the chiral-hexatic phase the nematic order 
freezes [|| and we can rewrite ([j]) in the "London limit" : 

F bulk = J d 3 x^K A (V9 6 ) 2 -K A q Q n Q - V0 6 j, (2) 

where Ka = 72\tpe\ 2 . In the following we take the average 
nematic axis to be no = z. This free energy, however, 
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applies only deep inside the sample. In general there 
will be surface terms which can bias the hexatic bond- 
order direction. For instance, if the boundary is a flat 
surface, then the mesogens will either be "surface loving" 
or "surface hating" and will try to maximize or minimize 
their contact with the surface, respectively. If the surface 
tension is large, the first few layers of molecules will be 
locked and then simple geometric arguments, akin to the 
classic Onsager treatment of nematics M , show that there 
is an entropic penalty for the hexatic order to rotate away 
from the preferred direction, as shown in Figure 0. 
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FIG. 2. Hexagonally close-packed disks near a wall with 
surface normal N. Each hexagonal crystallite defines a set 
of hexatic axes shown as six-pointed arrows. In (a) we show 
a favorable packing for "surface loving" disks, while in (b) 
we show a packing of "surface hating" disks. For either case, 
deviations away from this packing are disfavored entropically. 

To model the effect of this surface interaction, we add 
to (0) the surface term: 
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dsdz h cos(6#e) 
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where fi is the volume of the sample and dSl is the bound- 
ary in the xy-plane. The sign of h is determined by the 
surface affinity of the mesogens, while the magnitude of h 
can be estimated through an excluded- volume, entropic 
analysis. Taking the mesogens to be polymers (so that in 
Figure |[ the disks are the polymer cross-sections in the 
xy-plane), we consider an hexagonal crystallite of extent 
£ in the xy-plane and length Lp along z, rotating away 
from its preferred direction. Presumably Lp is precisely 
the polymer persistence length, which is the typical cor- 
relation of the polymer positions along the nematic axis. 
In this case, the entropy cost for a rotation 9 is 
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where 7 < 1 is a factor of order unity, N = V/(Lp£ 2 ) is 
the number of polymer crystallites and V is the volume 
of the system. Converting (Q) into an integral over the 



surface gives \h\ « k B T/(Lpt;). For the polymer hexatic 
phase studied in ||, L P m 500 A and £ w 400l. 

We propose that this surface anchoring can prevent the 
equilibrium twisting that a chiral-hcxatic phase must ex- 
hibit whenever the nematic director is aligned. To pur- 
sue this, we consider configurations in which the hexatic 
bond-order is uniform at any constant height z. This 
should be reasonable deep inside the hexatic phase where 
there is true long-range bond order. The free energy is: 
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where K = KaA, H = hL, A is the cross-sectional area 
of the sample and L is the length of the sample along the 
wall which provides the anchoring. This effective free 
energy has frequently been studied before in the con- 
text of commensurate-incommensurate transitions. For 
small values of go the cos(6#6) term prevents any twist- 
ing, while for larger values the system admits soliton so- 
lutions which let 6q slip by 2tt/6. In equilibrium, it is 
straightforward Go] to show that a soliton is the lowest 
energy configuration when 
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In addition to the estimate for h, we can estimate Ka 
via dimensional analysis: Ka ~ k B T/^. How should we 
estimate q ? In the cholesteric phase of DNA, the pitch 
is on the order of 1/j.m. We can use this to estimate the 
rate at which a polymer bundle twists about the aver- 
age nematic axis. To do so we consider a bundle of N 
polymers, with separation a. We let them twist about 
the center of the bundle so that the polymer that crosses 
z = at (p, 4>) has conformation [pTj 



H(z;p, 4>) = [pcos(q Q z + <fi),psm(q z 
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We calculate the energy for a twisted bundle of radius po 
and length L z . There are two contributions to the energy: 
the bending energy of the individual polymers and the 
Frank free energy of the liquid crystal deformation. The 
first energy is proportional to the bending modulus n. If 
the polymers sit at radii pj and have an average spacing 
a, then, for qop <C 1 we have: 
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where TV is the number of polymers in the bundle. The 
liquid crystal deformation has no splay and we have just 
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accounted for the bend deformations. The remaining en- 
ergetic contribution is from the twist: 



(2<7o + go) 



Finally, minimizing the total free energy F 
-Ftwist over q , we find: 
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The bending modulus is k = k B TLp where Lp is the 
persistence length, while K 2 — k B T/d where d is the 
Odijk length, given by d ~ (Lpa 2 ) 1 / 3 . Thus an esti- 
mate of the coefficient of eft is Nn 2 L 4 J 3 a 2 / 3 / P 2 where P 
is the hexatic pitch. The experiment || finds a trans- 
lational correlation length on the order of ten a = Anm 
intermolecular separations so presumably each twisting 
bundle is composed of N ~ 100 molecules. In this 
case, Nk/K2 ~ 10 _2 //m 2 . If we take a typical range 
1/im < 27r/<7o < 10/im then ( |l0| ) shows that 2//m < 
2ir/qo < 20/im. This follows from the fact that though 
the persistence length is long for DNA, it is still two 
orders of magnitude smaller than the cholesteric pitch: 
twisting on this length scale is easy for the polymers. . 

Using (H) , we find that twisting (through the introduc- 
tion of solitons in the 0q field) becomes unfavorable for 
systems of linear dimension smaller than a certain criti- 
cal length. For system sizes such that L < 16/ [LpTT 2 q 2 ^ 
the boundary terms dominate the free energy so the bond 
order is effectively pinned by these surface effects. From 
the above estimates of the equilibrium pitch of the hex- 
atic bond order we find that the largest system which 
allows surface pinning is on the order of 300/im when the 
pitch is 20/im. The sample was on the order of 1mm fl^ j 
so the outer boundaries of the sample are possibly too far 
away to cause lock-in. However, it is possible that sam- 
ple preparation introduced interior walls. A wall spacing 
on the order of ~ 300/im would imply internal regions 
of roughly 10 7 DNA molecules. If, on the other hand, 
we were to take the smallest reasonable pitch of 1/im, we 
would find a wall spacing of 3/im, which is unreasonable. 
Indeed, it is likely that in the hexatic phase, intermolec- 
ular correlations can modify the chiral interaction as ar- 
gued in H) and |l3j and lead to a lengthening of the pitch. 
It is thus reasonable that in the hexatic phase the chiral 
strength can be reduced and we are somewhere between 
these two extremes. 

The data, if taken to represent only the bulk physics, 
shows that 2n/qo ~ 1mm since the hexatic order does 
not twist by 2tt/6 inside the illuminated region. How- 
ever, with walls one millimeter away, the data is con- 
sistent with an equilibrium pitch of 44/im, not that far 
from our naive estimate of the pitch. Thus, bearing in 
mind the Landau free energy of this system, one can only 
conclude that the hexatic pitch is longer than 44/im, but 



not infinite! The absence of the observation of a finite 
hexatic pitch can be attributed to the surface pinning 
effects. If we take our estimate of the hexatic pitch and 
make the reasonable assumption of a few internal walls 
in the sample, or if we imagine the hexatic pitch to be 
within a factor of two of typical pitches, it is clear that 
such a pinned state is expected. In either case the effects 
of surface pinning are required to interpret the data. As 
further evidence for our pinning hypothesis, it is known 
that the hexatic director tends to align with the long axis 
of the sample [fL2| . 




-1 1 
z (|im) 



FIG. 3. Three solutions for the bond-order field as a func- 
tion of height. We show the one soliton {k = 1, the dotted 
line) solution, as well as multi-soliton solutions for k < 1. 
The parameter k is determined by the sample dimension L 
and the equilibrium pitch go allowing the construction of the 
phase diagram shown in the inset. In the "pinned" phase (P) 
the bond-order field 9e is prevented from twisting by the sur- 
face pinning while in the "free" (F) phase the bond-order can 
twist. Upon approaching the phase boundary (solid line) from 
below, the bond order begins to twist by the creation of soli- 
tons. The actual phase boundary line is chosen by arbitrarily 
assigning a maximum (non-zero) soliton density as described 
in the text. 

Since both sides of (||) are the same order of magnitude 
the condition to prevent solitons is, at best, marginally 
satisfied. What happens if we have solitons? We shall 
see that even close to the transition a state with solitons 
will still have an achiral hexatic signature. When (^) 
is an equality, it becomes energetically favorable for the 
system to introduce one soliton into the bond order field, 
0q(z). Such a field configuration, in which the hexatic 
bond order jumps by 27r/6 near z = 0, is shown by the 
dotted line in Figure ||. At low densities two solitons 
interact via a weak exponential repulsion with a range 
set by the distance over which 6q (z) changes rapidly with 
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This weak interaction of two solitons leads to the rapid 
proliferation of such twists in the system upon tuning the 
system parameters, L and qo. Once these parameters 
are tuned so that the chemical potential of a soliton be- 
comes negative, numerous solitons are spontaneously cre- 
ated until their weak mutual repulsion returns the soliton 
chemical potential to zero. In this case the bond-order 
field takes the form ifl: 
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where am(i, to) is the inverse of the incomplete elliptic 
integral of the second kind and the parameter k is related 
to the energy via E(k 2 ) — kirq^^J LpL/ A where E(x) is 
the complete elliptic integral of the second kind. The 
solution to ( |ll|) describes a regular array of solitons with 
a density controlled by the value of the parameter k. 

As shown in Figure |] the parameter k controls the equi- 
librium density of solitons in the system. When k = 1, 
there is precisely one soliton in the infinite system and, as 
k decreases from unity, the soliton density grows mono- 
tonically. When k — 0.95 this soliton density reaches a 
point where mutually overlapping solitons join to form 
an almost uniform gradient in the bond order. The sys- 
tem returns to its "free" behavior with 6q (z) oc z. On 
the other hand, if k is larger than unity, there are no soli- 
tons in the system in equilibrium. The bond-order field 
is completely pinned by the edges of the sample so that 
9q (z) is constant. 

The experimental signature of the pinned phase is a 
six-fold modulation of the in-plane scattering intensity. 
It can, in fact, persist into regions of the phase diagram 
which allow modest densities of solitons. The six- fold 
modulation of the scattering intensity should remain ev- 
ident as long as the density of solitons is small enough so 
that distinguishable plateaus in the 9q shown in Figure 
H exist. We take a somewhat arbitrary but nevertheless 
conservative estimate of the transition to be the point 
at which the plateaus are twice as long as the interven- 
ing regions where the Q§ field changes rapidly |l5| . By 
this criterion the solid curve in Figure [| represents a free 
phase configuration of 6§ . The division of the parameter 
space between the pinned and freely-twisting bond order 
selects a critical value of k = k* . Solving for k* leads 
to a phase boundary between pinned and free states of 
the form: qoL 1 / 2 = const. This curve (with the constant 
set by k*) is shown in the inset of Figure ||. Addition- 
ally, the soliton lattice could, in principle, be observed 
via scattering due to its periodic structure along z. 

Finally, we comment on the role of the thermal effects. 
Fluctuation effects become important near the boundary 
of the pinned and free phases since there the system ex- 
hibits a delicate balance of the bulk elastic and surface 
energies. The qualitative effect of thermally produced 



solitons is to move the phase boundary down and to the 
left in the inset of Figure |L To assess the magnitude 
of this shift, we compare the surface free energy cost to 
kftT '. In a system with a size as small as L ~ 10/im the 
surface energy is on the order of 10 3 /cbT- The transi- 
tion is controlled then by the balance of energies three 
orders of magnitude greater than the thermal energy, so 
we conclude that thermal effects are negligible except in 
the narrow region along the phase boundary where the 
surface and bulk energies are balanced to within 0.1%. 
Additionally, at higher temperatures the long-range hex- 
atic order will melt, leaving a surface hexatic along with 
a bulk cholesteric. The details of this transition could be 
confirmed through further experiments along the lines of 
1- 
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